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ABSTRACT 

Rich and massive clusters of galaxies at intermediate redshift are capable of magni¬ 
fying and distorting the images of background galaxies. A comparison of different mass 
estimators among these clusters can provide useful information about the distribution 
and composition of cluster matter and their dynamical evolution. Using a hitherto 
largest sample of lensing clusters drawn from literature, we compare the gravitat¬ 
ing masses of clusters derived from the strong/weak gravitational lensing phenomena, 
from the X-ray measurements based on the assumption of hydrostatic equilibrium, 
and from the conventional isothermal sphere model for the dark matter profile char¬ 
acterized by the velocity dispersion and core radius of galaxy distributions in clusters. 
While there is an excellent agreement between the weak lensing, X-ray and isothermal 
sphere model determined cluster masses, these methods are likely to underestimate 
the gravitating masses enclosed within the central cores of clusters by a factor of 2-4 as 
compared with the strong lensing results. Such a mass discrepancy has probably arisen 
from the inappropriate applications of the weak lensing technique and the hydrostatic 
equilibrium hypothesis to the central regions of clusters as well as an unreasonably 
large core radius for both luminous and dark matter profiles. Nevertheless, it is pointed 
out that these cluster mass estimators may be safely applied on scales greater than the 
core sizes. Namely, the overall clusters of galaxies at intermediate redshift can still be 
regarded as the dynamically relaxed systems, in which the velocity dispersion of galax¬ 
ies and the temperature of X-ray emitting gas are good indicators of the underlying 
gravitational potentials of clusters. 

Key words: dark matter - galaxies: clusters: general - gravitational lensing - X-rays: 
galaxies 


1 INTRODUCTION 

Clusters of galaxies are the largest coherent and gravita¬ 
tionally bound objects in the universe. They are often used 
for cosmological test of theories of structure formation. In 
particular, their matter composition, the mass-to-light ratio 
M/L (e.g. Bahcall, Lubin, & Dorman 1995) and the baryon 
fraction ft, (e.g. White et al, 1993) play a potentially impor¬ 
tant role in the direct measurement of mean mass density of 
the universe, U m . It appears that the observational and the¬ 
oretical results are likely to merge nowadays, in favor of a low 
mass density universe of ~ 0.3 (e.g. Bahcall et al. 1995; 
Ostriker & Steinhardt 1995). Yet, the key issue in such a “di¬ 
rect” measurement of the cosmological density parameter is 
closely connected to the question of how accurately one can 
determine the total gravitating masses of galaxy clusters. 


Recall that the traditional cluster estimators strongly rely 
upon the assumption of hydrostatic equilibrium: both opti¬ 
cal galaxies and intracluster diffuse gas quasi-statically trace 
the underlying gravitational potential of the whole cluster. 
However, a number of recent X-ray observations (e.g. Henry 
& Briel 1995; Markevitch 1996; etc.) have argued that the 
hydrostatic hypothesis is inapplicable to at least merging 
clusters. Therefore, an independent method of estimating 
cluster mass is highly appreciated. At present, it is widely 
recognized that the gravitational lensing is a unique and also 
ideal tool to fulfill the task, in the sense that it provides a fair 
estimate of cluster masses regardless of the cluster matter 
components and their dynamical states. 

In the framework of gravitational lensing, mapping the 
gravitational potential of a galaxy cluster can be achieved 
by the detailed modeling of the strongly distorted and mag- 
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nified arclike images inside the core of the cluster (Fort & 
Mellier 1994), by the statistical analysis of the weakly dis¬ 
torted images of faint and distant galaxies around the clus¬ 
ter (Kaiser & Squires 1993), and by counting the popula¬ 
tions of distant galaxies (Wu & Hammer 1995; Taylor et al. 
1998) and quasars (Wu & Fang 1996a; references therein) 
around the cluster. A combination of these lensing phenom¬ 
ena would, in principle, allow one to depict an overall matter 
distribution of cluster on scales spanning two decades from 
the inner core of ~ 30 kpc to the outer regime of ~ lOMpc 
(Wu & Fang 1996b). However, in the practical measurements 
the lensing signals often appear to be relatively weak, and a 
precise determination of cluster masses suffers from the limi¬ 
tations of instrument sensitivities and viewing fields. In par¬ 
ticular, for different lensing phenomena different techniques 
have been developed and employed in the reconstruction of 
matter distributions of clusters. Hence, before one proceeds 
to the cosmological applications of the gravitating masses of 
clusters derived from gravitational lensing, it is also worthy 
of examining whether different approaches based on differ¬ 
ent lensing phenomena yield a consistent gravitating mass 
of cluster. 

On the other hand, apart from the doubt as to whether 
the gravitational lensing provides an accurate cluster mass 
estimate, it is of great interest to compare the cluster masses 
derived from traditional method with those from gravita¬ 
tional lensing. Previous comparison has been made among 
many individual clusters, in which both strongly/weakly 
distorted images of background galaxies and X-ray/optical 
emissions are detected. Basically, a good agreement between 
the X-ray and weak lensing determined cluster masses has 
been detected (Squires et al. 1996, 1997b; Smail et al. 1997; 
Allen 1998), in contrast to the situation of strong lensing, in 
which the lensing derived cluster masses are typically ~ 2 —4 
times larger than the cluster dynamical masses obtained un¬ 
der the assumption of hydrostatic equilibrium (e.g. Wu 1994; 
Miralda-Escude & Babul 1995; Allen, Fabian, & Kneib 1996; 
Markevitch 1997). This argument is further strengthened by 
Wu & Fang (1997) using a large sample of 29 lensing clusters 
drawn from literature. Several attempts have correspond¬ 
ingly been made to explain the discrepancy such as the possi¬ 
ble contributions of nonthermal pressure (Loeb & Mao 1994; 
Ensslin et al. 1997) and projection effect (Miralda-Escude & 
Babul 1995; Girardi et al. 1997b). An important progress in 
understanding the issue has been made recently by Allen 
(1998), based on a detailed comparison of the cluster mass 
determinations from X-ray and gravitational lensing for 13 
lensing clusters. He showed that the mass discrepancy be¬ 
tween the X-ray and strong lensing analyses exists only in 
the non-cooling flow clusters, which can be well accounted 
for by the significant offsets between the X-ray and strong 
lensing centers. This implies that in the central regions of 
these clusters, the traditional “quasi-static” hypothesis for 
the X-ray emitting gas has probably broken down. However, 
it should be noticed that even for these non-cooling clusters 
the weak lensing determined cluster masses agree nicely with 
the X-ray determined ones. 

The current comparisons of different cluster mass esti¬ 
mates seem to suggest that the reported mass discrepancy 
between X-ray and gravitational lensing analyses may, at 
least partially, arise from the employments of different lens¬ 
ing techniques: Strong lensing and weak lensing yield an in¬ 


consistent cluster mass. Although this disagreement could be 
attributed to observations because strong lensing and weak 
lensing measurements probe the different regions of clus¬ 
ters, a direct comparison between cluster masses obtained 
from these two lensing phenomena has now become possible 
with the progress of weak lensing technique and the increas¬ 
ing population of arcs/arclets. If the cluster masses derived 
from strong lensing do systematically exceed the weak lens¬ 
ing values, one would be faced with the following difficulties: 
The consistency of the weak lensing and X-ray determined 
cluster masses indicates that the strong lensing analysis may 
lead to an overestimate of gravitating mass for cluster, which 
requires a re-examination of the modeling of arcs/arclets as 
the cluster mass estimator. On the other hand, the strong 
lensing phenomena (arcs or arclets) are usually believed to 
be a good indicator of the underlying gravitational poten¬ 
tial of the lensing clusters, in which there are almost no 
free parameters once the redshifts of the arclike images are 
known. This leads us to the possibility that both weak lens¬ 
ing and X-ray measurements may underestimate the true 
cluster mass by a factor of 2-4. The latter point is directly 
related to the question whether the previously estimated 
M/L and fb of clusters gives rise to an underestimate of the 
cosmological density parameter fl m . Of cause, there may ex¬ 
ist a third possibility that enables us to reconcile the cluster 
mass difference while maintain the validity of these differ¬ 
ent mass estimators. Alternatively, an examination of the 
consistency or discrepancy between the strong/weak lensing 
and X-ray/optical mass measurements of clusters will be of 
great help to our understanding of dynamical and evolu¬ 
tionary properties of rich clusters at intermediate redshifts. 
Throughout this paper, we assume Ho = 50 km s” 1 Mpc -1 
and flo = 1. 


2 SAMPLE 

By searching the literature, we have collected a strong lens¬ 
ing cluster sample of 38 clusters (Table 1), which contains 48 
arc/arclet images of distant galaxies. We adopt only the pub¬ 
lished data sets of the projected cluster masses within the 
positions of arcs/arclets and make no attempt as far as pos¬ 
sible to extrapolate the original work. Therefore, a number 
of fainter and smaller arclets have not been included in Table 
1. Our weak lensing cluster sample consists of 24 clusters, 
in which the projected radial mass distributions have been 
given for 10 clusters while only one mass measurement inte¬ 
rior to a certain radius is listed for the remaining clusters. 
It is helpful to recall how the cluster mass is computed in 
the framework of gravitational lensing. In the case of strong 
lensing, the projected cluster mass within the position of 
arc/arclet r ar . c is often estimated through 

1Tllens,arc (v <C Tare) — T^T arc ^j C rit , (1) 

where £ C rit = (c 2 / A-kG)(D s / DdDds) is the critical surface 
mass density, with Dd , D s and Dds being the angular di¬ 
ameter distances to the cluster, to the background galaxy, 
and from the cluster to the galaxy, respectively. Eq.(l) is 
actually the lensing equation for a cluster lens of spherical 
mass distribution with a negligible small alignment param¬ 
eter for the distant galaxy as compared with r a rc- As for 
weak lensing, a similar expression to eq.(l) is applied: 
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Table 1 . Strong Lensing Cluster Sample 


cluster 

Zcluster 

T (keV) 

<?gal (km S J ) 

Zarc 

rar-c(Mpc) 

?nass(lO 14 M0) 

A370 

0.373 

7 1+ 1 - 0 

13671?“ 

1.3 

0.35 

13.0 





0.724 

0.16 

2.90 

A963 

0.206 

ft q+ 0.4 
b -^-0.5 

nooj^g 


0.0517 

0.25 





0.711 

0.080 

0.60 

A1204 

0.170 

q 0 + 0.2 
^ — 0.2 

19891 245 


0.0311 

0.10 

A1689 

0.181 

9 0 " 1 ” 0 ' 4 
y * u -0.3 


0.183 

3.6 

A1835 

0.252 

8 . 2 t °; 5 5 



0.150 

1.98 

A1942 

0.224 




0.0372 

0.18 

A2104 

0.155 




0.025 

0.064 

A2163 

0.203 

i3-9tg:S 

1680 

0.728 

0.0661 

0.41 

A2218 

0.171 

7.0±i;° 

14051“? 

1.034 

0.26 

2.7 





0.702 

0.0794 

0.623 





2.515 

0.0848 

0.570 

A2219 

0.228 

n- 8 ±i;l 



0.079 

0.517 






0.110 

1.60 

A2280 

0.326 

8 9 +1 ' 0 
°- y - 0.8 

948l 5 2 “ 


0.080 

0.59 

A2390 

0.228 

10931?? 

0.913 

0.177 

2.54 

A2397 

0.224 




0.0698 

0.45 

A2744 

0.308 

12 1+ 1 - 4 

19501??? 


0.1196 

1.136 

A3408 

0.0419 



0.0728 

0.0568 

0.44 

S295 

0.299 


907 


0.0329 

0.14 





0.930 

0.14 

1.6 

CL0024 

0.391 


1339+ 233 
xoo»_ 233 

1.390 

0.214 

3.324 

CL0302 

0.423 


1100 


0.122 

1.60 

CL0500 

0.327 

7 9+3.7 

11521??? 


0.15 

1.90 

CL2236 

0.552 



1.116 

0.0876 

0.30 

CL2244 

0.328 

6 >i:l 

6061 ?? 

2.236 

0.0465 

0.20 

MS0440 

0.197 

5-3 tli 

0.530 

0.089 

0.89 

MS0451 

0.539 

io-2±i:| 

13711??? 


0.190 

5.2 

MS0955 

0.145 




0.0385 

0.22 

MS1006 

0.261 


9061??? 


0.079 

0.57 






0.14 

1.8 






0.28 

7.2 

MS1008 

0.306 

7 o+2.5 

1054+ 107 
±uo — 107 


0.26 

6.1 

MS1137 

0.783 




0.044 

0.19 






0.147 

2.1 






0.151 

2.2 

MS1358 

0.329 

ft r+0.7 
d, 0 -0.6 

9371?? 

4.92 

0.121 

0.827 

MS1445 

0.257 

D,D -0.3 

H33 + 140 

HOO_140 


0.098 

0.86 

MS1621 

0.427 


7931?? 


0.046 

0.24 

MS1910 

0.246 




0.33 

9.6 

MS2053 

0.523 




0.119 

2.60 

MS2137 

0.313 

4 4+0-4 
^-0.4 

960 


0.0874 

0.71 

MS2318 

0.130 

5.1 



0.12 

1.3 

AC 114 

0.310 

8-+0.9 

16491??? 

0.639 

0.35 

13.0 

GH02154 

0.320 



0.721 

0.0342 

0.20 

PKS0745 

0.103 

8-5ti1 


0.433 

0.0459 

0.30 

RXJ 1347 

0.451 

ii 

1235 


0.24 

4.2 


Data are collected from Wu & Hammer (1993), Le Fevre et al. (1994), Kneib Sz Soucail (1995), Wu & Fang (1997), Allen (1998), 
Campusano, Kneib & Hardy (1998), Tyson, Kochanski & Dell’Antonio (1998) and Clowe et al. (1998) 
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Table 2. Weak Lensing Cluster Sample 


cluster 

% cluster 

T (keV) 

<?gal (km S X ) 

r(Mpc) 

mass(10 14 MQ) 

references 

A1689 

0.1810 

9-0tS1 

1989^245 

0 . 20 - 1.6 

3.2-14 

Tyson et al. (1995) 

A2163 

0.2030 

i3-9±S:S 

1680 

0.09-0.89 

0.34-4.05 

Squires et al. (1997b) 

A2218 

0.1710 

7 0 + 1 ' 0 

1 - u -i.o 

1405 +163 

iiuo_i4 5 

0.8 


Kneib et al. (1995) 





0.405 

2 -lOt 3 ; 3 ® 

Smail et al. (1997) 





0.25-1.16 

1.37-11.7 

Squires et al. (1996a) 

A2263 

0.208 



0.87 

5-o ± 2 :| 

Allen (1998) 

A2390 

0.2279 

8 9 + 1 ' 0 
°- a - 0.8 

1093+® 4 

0.19-1.20 

1.59-19.1 

Squires et al. (1996) 





0.94 

toll 

Allen (1998) 

A2744 

0.308 

12 1+ 1 ' 4 
- 1.0 

1950^334 

0.40 

3 70 +°' 64 
(u _0.64 

Smail et al. (1997) 

0957+561 

0.36 


71 c + 130 

1 lo -130 

0.06-0.80 

0.12- 1.5 

Fischer et al. (1997a) 

3C324 

1.206 



0.5 

6.0 

Smail et al. (1995) 

3C295 

0.46 

12.6 

1670”*” 364 

io/u_ 364 

0.40 

4 70+°- 76 
4 ' ,U -0.76 

Smail et al. (1997) 

3C336 




0.50 

4 8+ 1 ' 0 

Bower et al. (1997) 

AC114 

0.3100 


1649+^0 

0.50 

4 0 +°' 4 

^■ u -0.4 

Allen (1998) 

CL0016 

0.5545 

8 . 0 ±i;° 

1234+ 4 28 

0.40 

3 74+ 1 - 2 ® 

Smail et al. (1997) 

CL0024 

0.3910 


1339+ 233 
i0jy _ 233 

3.0 

40 

Bonnet et al. (1994) 




0.40 

2 7fi +°- 74 

1 -0.74 

Smail et al. (1997) 

CL0054 

0.56 



0.40 

3 42+ 1 - 2 ® 

Smail et al. (1997) 

CL0303 

0.0349 


1079j:235 

0.40 

n 44 + 0 - 90 
u -^-0.90 

Smail et al. (1997) 

CL0412 

0.51 



0.40 

°- 5 °ia50 

Smail et al. (1997) 

CL0939 

0.4510 


lOSlli ® 4 

0.75 

6-i 

Seitz et al. (1996) 





0.40 


Smail et al. (1997) 

CL1601 

0.54 


1166 

0.40 

1 -5 4 ±i;|| 

Smail et al. (1997) 

MS1054 

0.826 

14-7 ttt 

1643t®06 

0.23-2.0 

0.62-27.3 

Luppino et al. (1997) 

MS1137 

0.783 



0.18-1.21 

1.60-4.74 

Clowe et al. (1998) 

MS1224 

0.3255 

43 + 115 

802 ±1° 

0.96 

7.0 

Fahlman et al. (1994) 

MS1358 

0.3290 

r+0.7 

D *°-0.6 

937+g 4 

0.12-1.29 

0.38-3.65 

Hoekstra et al. (1998) 

RXJ1347 

0.4510 

11 4 +1,1 

1235 

0.24-2.6 

2.6-30 

Fischer et al. (1997b) 

RXJ1716 

0.813 

6 7 + 2,0 
* — 2.0 

1892 

0.18-1.11 

0.95-5.90 

Clowe et al. (1998) 


177'lens, weak (r < ro) = 7rroC(ro)£cri*. (2) 

Here £ is defined as (Fahlman et al. 1994) 

CM = 2 (l - ro/r-jL*) 1 fr™° X ('y T )d In r 

= ff(r < r 0 ) - a(r 0 < r < r max ), (3) 

in which {yr) is the mean tangential shear around a circular 
path of radius r introduced by the gravitational potential of 
the cluster, a(< r o) and <j(ro < r < rmax) represent the 
mean surface mass density interior to ro and in the annulus 
ro < r < rmax, respectively. Therefore, eq.(2) provides a low 
bound on the projected cluster mass within ro. 

In order to reconstruct the cluster mass rni erlSiarc (r < 
r a rc ) or miens,weak(r < ro), one needs also to know the red- 
shift of the arc/arclet or the spatial distribution of back¬ 
ground galaxies. Spectroscopic measurements have been 
made for about ~ 2/5 of the arcs/arclets in Table 1. When 
the data of redshifts are not available, a redshift of z s = 0.8 
is often used in the estimate of E cr .jt, except for the arcs in 
clusters MS1137+66 at z = 0.783 where z s = 1.5 is adopted 
in Table 1. Such an assumption is justified by the recent 


work of Ebbels et al. (1998), who obtained the spectroscopic 
identifications of 18 arclike images behind A2218 and found 
a mean redshift of (z a ) = 0.8-1 at R ~ 25.5. Nevertheless, 
placing the background galaxies at z s = 2 instead of z s = 0.8 
would reduce the value of E cr .jt by a factor of only 1.4 for a 
mean cluster redshift of (zd) ~ 0.3. Weak lensing measure¬ 
ments involve a great number of faint and distant galaxies 
behind foreground clusters, and the spectroscopic measure¬ 
ments of all the population of galaxies in the fields seem 
to be impossible at present. Different calibrations have thus 
been applied by different authors in the determinations of 
Ecr-it- For clusters at relatively lower redshifts of Zd ~ 0.3, it 
appears to be plausible to assume the mean redshift of the 
background galaxies to be at z s ~0.8-1.2 when combined 
with the existing work of deep galaxy surveys. However, it 
should be noted that the big uncertainty of up to a factor 
of ~ 5 arising from the unknown redshift distributions of 
background galaxies may occur in the evaluation of E cr n 
for high-3 clusters at Zd « 0.8 (e.g. Luppino & Kaiser 1997). 

Before we turn to the comparison of strong lensing, 
weak lensing and X-ray/optical determined cluster masses, 
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we display in Fig.l the relationship between the velocity 
dispersion (a ga i) of galaxies and the temperature (T) of X- 
ray emitting gas for 20 lensing clusters in Table 1 and Ta¬ 
ble 2, for which both a ga i and T are available in literature. 
It is believed that the o ga i — T relationship can provide a 
straightforward yet robust test for the dynamical properties 
of clusters of galaxies (Cavaliere & Fusco-Femiano 1976). 
Employing the least-square fit of a power-law to the data of 
Fig.l yields (ag„xl km s" 1 ) = 10 2 - 57±013 (fcT/keV)°' 59±014 . 
Regardless of its large error bars, this relationship is essen¬ 
tially identical with the one for the nearby clusters (Girardi 
et al. 1996 and references therein), and is also consistent with 
the isothermal hydrostatic scenario of a ga i ~ T°’ 5 . Based on 
the o' g ai — T relationship alone, we may conclude that no sig¬ 
nificant dynamical and cosmic evolution has been detected 
for those massive lensing clusters within intermediate red- 
shift (zd) ~ 0.3. It is interesting to recall that similar con¬ 
clusions have been reached by numerous recent studies on 
X-ray clusters. For instances, the cluster number counts ex¬ 
hibit no evolutionary tendency at least to redshift of as high 
as z ~ 0.8 (Rosati et al. 1998), and no significant differences 
in the X-ray luminosity-temperature relationship and the 
velocity dispersion-temperature relationship between low- 
redshift and high-redsliift clusters are seen (e.g. Mushotzky 
& Scharf 1997). In particular, the distribution of core ra¬ 
dius of the intracluster gas of nearby clusters accords with 
that of distant clusters (z > 0.4) (Vikhlinin et al. 1998). So, 
these arguments may eventually support the hypothesis that 
overall, both galaxies and gas are the tracers of the depth 
and shape of the underlying gravitational potential of clus¬ 
ter, in despite of the presence of substructures and merging 
activities on small scales. 


3 COMPARISONS OF DIFFERENT MASS 
ESTIMATES 

3.1 Strong and weak lensing 

Ten weak lensing clusters in Table 2 also contain arclike im¬ 
ages of background galaxies in their central cores. Compar¬ 
ison of strong and weak lensing determined cluster masses 
among these clusters are straightforward and shown in Fig.2. 
Except for A1689 where the absolute mass calibration was 
made using the giant arcs which trace the Einstein radius 
(Tyson & Fischer 1995), the rest strong lensing events seem 
to yield larger cluster masses than the weak lensing mea¬ 
surements when the two methods become to be comparable 
with each other at the central regimes. In Fig.3 we display 
the strong/weak lensing determined cluster masses at differ¬ 
ent radii utilizing all the data sets in Table 1 and Table 2. 
Although such a comparison sounds less serious in the sense 
that they are different clusters, it is by no means of physical 
insignificance. Indeed, the majority of clusters of galaxies 
that are capable of magnifying and distorting the images of 
background galaxies are very rich clusters at intermediate 
redshift. They are representative of the (most) massive clus¬ 
ters at (z) ~ 0.3. Statistically, their gross dynamical proper¬ 
ties should look very similar. It is remarkable that the data 
points in Fig.3 are clearly separated into two parties: the 
strong lensing determined cluster masses systematically and 
significantly exceed the weak lensing values. Actually, the 


strong lensing data appear in the plot as if they were the 
upper limits to the weak lensing measurements. Alterna¬ 
tively, the cluster masses given by the arcs/arclets with and 
without confirmed redshifts show no significant differences. 
The very small dispersion in the strong lensing results arises 
from the fact that the critical surface mass density parame¬ 
ter Tj cr u remains roughly unchanged for the known lensing 
systems which have a mean cluster redshift of (zd) ~ 0.3 
and a mean source redshift of (z s ) « 0.8. Consequently, the 
“mass density profile” derived from the strong lensing events 
varies as r~ 2 . 


3.2 Lensing and X-ray measurements 


X-ray observations have been made for most of the lensing 
clusters. Under the assumption that the hot diffuse gas is 
in hydrostatic equilibrium with the underlying gravitational 
potentials of clusters, one can easily obtain the X-ray cluster 
masses provided that the gas and temperature radial profiles 
are well determined. Adopting the conventional isothermal 
P model for gas distribution which is characterized by the 
core radius r xra y,c, the index /3fu and the temperature T, 
we can write out the projected X-ray cluster mass within 
radius r to be (Wu 1994) 


m X ray = 1.13 x 10 13 /3/um(r) 


where 
m{r) = 




R/ r xray ,c 
ray ,c 


1 f'xrayjC \ 

( kT \ 

yo.i Mpc J 

U keVj 

,e) 3 





r / r xr 


r \2 3+ar 

• xray,c) (i _|- x ^)^ a " X/ ’ 


and R is the physical radius of the cluster and will be taken 
to be R = 3 Mpc in the actual computation. Our conclusion 
is unaffected by this choice. 

We have computed the projected X-ray cluster mass 
m X ray interior to the position of each arclike image r arc 
or the corresponding radius of each weak lensing measure¬ 
ment r o for the 21/13 strong/weak lensing clusters of known 
temperatures listed in Table 1 and Table 2, and the result¬ 
ing rrixray versus mi e „s are plotted in Fig.4 for fifu = 2/3 
and a mean core radius of ( r X ra y ,c ) = 0.25 Mpc. It turns 
out that although both the strong and weak lensing de¬ 
termined cluster masses miens,arc and miens,weak show a 
good correlation with the X-ray masses m xray , their am¬ 
plitudes are very different: miens,arc/m xra y = 3.23 ± 1.21 
and Trilens,weak / m xra y = 0.97 ± 0.44. The scatter in the fit 
of miens,weak /m xra y can be reduced when cluster A2163 is 
excluded (Fig.4b) (Note that A2163 is one of the hottest 
clusters known so far). That is to say, there is a signifi¬ 
cant discrepancy between the strong lensing derived clus¬ 
ter masses and the X-ray cluster masses, while an excellent 
agreement between the weak lensing results and the X-ray 
masses is found. The only way to reconcile miens,arc with 
m xra y is to adopt a considerably small core radius for the 
gas profile. Fig.5 demonstrates another plot of miens against 
m xr ay by reducing r xra y,c to 0.025 Mpc. In this case, the ra¬ 
tios of miens,arc and miens,weak to m xr a y read 1.42 ± 0.87 
and 0.91 ± 0.35, respectively. Yet, this does not relax the 
disagreement between miens,arc and miens,weak- Rather, the 
employment of a smaller core radius leads to a significant 
increase of the X-ray cluster mass estimate at small radius. 
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Recall that the crucial point behind the remarkable agree¬ 
ment of the strong lensing and X-ray determined masses for 
the cooling flow clusters reported by Allen (1998) is the small 
core radii of r X ray,c ~ 40-75 kpc, in contrast to r xray ,c « 250 
kpc for the non-cooling flow clusters. 

3.3 Velocity dispersion of galaxies as the tracer of 
cluster potential 

If galaxies trace the underlying gravitational potential of 
cluster, their velocity dispersion a ga i would be a good indi¬ 
cator of dark matter. A strict way to derive the dynamical 
mass of cluster from the distributions of galaxy population 
and their velocity dispersion is to work with the Jeans equa¬ 
tion. Here, we utilize an alternative approach to the issue. 
We adopt the so-called softened isothermal sphere model 
with a core radius r dar k,c for the total mass distribution of 
cluster, which is characterized by the velocity dispersion of 
galaxies. The original motivation was to examine whether 
(j gai provides a proper cluster mass estimate (Wu & Fang 
1997). The projected cluster mass within radius r is simply 

m op t(< r) = (^Jr 2 + r 2 arkc - r dark , c ) . (5) 

There are 21/18 clusters in Table 1/2 whose velocity disper¬ 
sions are observationally determined. We compute their to¬ 
tal masses in terms of eq.(5) at the corresponding positions 
of arclike images or weak lensing measurements. We first 
adopt a core radius of r dar k, c = 0.25 Mpc, in accord with the 
distributions of luminous matter (galaxies and gas) in clus¬ 
ter, and the resulting m op t are shown in Fig.6. It turns out 
that the weak lensing determined cluster masses mi en s,weak 
are in fairly good agreement with m op t despite of the large 
scatters: miens,weak/m op t = 1.08 ±0.70, whereas the strong 
lensing results miens, arc depart apparently from the expecta¬ 
tion of TTllens,arc ~ m,opt with TTllens,arc/TTlopt — 6.07 ± 3.98. 
Motivated by the argument that the dark matter profile is 
sharply peaked at the cluster center relative to the lumi¬ 
nous matter distribution (e.g. Hammer 1991; Wu & Ham¬ 
mer 1993; Durret et al. 1994; etc), we also present in Fig.6 
the results m op t for a smaller core radius of r dar k, c = 0.025 
Mpc. This indeed reduces the difference between mi en s,arc 
and m op t, but simultaneously breaks down the accordance 
of miens,weak and m 0 pt ■ The mean ratios of strong and weak 
lensing results to the isothermal sphere model determined 
cluster masses are now miens,arc/m op t = 1-44 ± 0.97 and 
miens,weak/rn op t = 0.63 ± 0.35, respectively. The scatter in 
miens,weak / m op t can be significantly reduced if the data sets 
of A2163 and RXJ1716 are excluded. 

A more reasonable way to estimate the dynamical mass 
of cluster in terms of velocity dispersion of galaxies is to 
employ the virial theorem: M v i r = 3a gal r v /G, which mea¬ 
sures the total cluster mass enclosed within a sphere of the 
so-called virial radius r v . A direct comparison of the virial 
and lensing cluster masses is somewhat difficult because the 
virial radius r v is usually much larger than the size which can 
be reached by the current gravitational lensing techniques. 
We display in Fig.7 the cluster masses given by strong/weak 
lensing and virial theorem for 10 clusters with the available 
data of Mnir in literature. The projected cluster mass de¬ 
rived from gravitational lensing would approach to the 3-D 
virial mass at radius as large as r v . However, the fact that 


the cluster regimes probed by the two methods show no over¬ 
laps can only lead us to arrive at the conclusion that these 
two independent methods seem to provide a consistent ra¬ 
dial matter distribution of cluster. 


4 DISCUSSION 

4.1 Strong lensing: overestimates cluster mass ? 

The consistence of the u ga i-T relationship of the lensing 
clusters with u ga i oc T 0 ' 5 expected under the hypothesis 
of isothermal and hydrostatic equilibrium suggests that the 
gas and galaxies are good tracers of underlying gravitational 
potential of the clusters. Indeed, a number of recent stud¬ 
ies have shown that clusters of galaxies suffer from little 
evolution, and their dynamical properties have remained al¬ 
most unchanged since 2 ~ 0.8 (Mushotzky, & Scharf 1977; 
Bahcall, Fan, & Cen 1997; Henry 1997; Rosati et al. 1998; 
Vikhlinin et al. 1988; etc.). This essentially justifies the em¬ 
ployment of the Jeans equation for X-ray emitting gas in 
the lensing clusters at intermediate redshifts. On the other 
hand, the weak lensing probes the gravitational potential 
fields of the clusters in a completely different way. So, the 
excellent agreement between weak lensing and X-ray deter¬ 
mined cluster masses serves as another convincing evidence 
for the lensing clusters being the dynamically relaxed sys¬ 
tems. 

A further examination of the mass determinations with 
eq.(l) tells us that there are no free parameters with which 
one can play when the redshifts of the lensing clusters and 
of the arclike images are observationally measured. It is un¬ 
likely that the replacement of the simple spherical lens model 
by a more sophisticated one can make a significant difference 
[see Allen (1998) for a detailed discussion; and references 
therein). Therefore, we are forced to accept the fact that 
the present modeling of strong lensing event does provide a 
reliable cluster mass estimate. 

However, from the comparisons of different cluster mass 
estimators in the above section, it is apparent that the clus¬ 
ter mass revealed by strong lensing exceeds that inferred 
from weak lensing and X-ray method by a factor of ~ 2 - 4. 
If we choose to trust the strong lensing result, we may need 
to break down the fairly good agreement between the weak 
lensing and X-ray derived cluster masses, and work out a 
mechanism which can account for the mass discrepancy be¬ 
tween the strong lensing and the latter two methods. There 
are quite a number of mechanisms such as projection effect 
and asymmetrical mass distributions which may result in an 
overestimate of cluster mass from arclike images. But, these 
mechanisms will simultaneously influence the weak lensing 
results. In the following discussion we will focus on whether 
weak lensing and X-ray analysis give rise to an underesti¬ 
mate of the gravitating masses of clusters. 

4.2 Weak lensing: underestimates cluster mass ? 

In principle, the weak lensing method eq.(2) and (3) always 
provides a lower bound on the cluster mass miens,weak (< ro) 
interior to radius ro- So, the total mass within ro could be 
considerably underestimated unless the outer control an¬ 
nulus r m ax is set to be sufficiently large and/or the true 
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mass density profile drops sharply along outward radius. In 
the current measurements of weak lensing effects, a value 
of Tmax ~ 1 — 2 Mpc is often adopted due to the limited 
data coverages. If cluster exhibits rather a flat surface mass 
density, e.g. a power-law of oc r 1 with index 7 < 1 , it is 
not impossible that the true cluster masses can be under¬ 
estimated by a factor of as large as ~ 2 , according to the 
formalism of Kaiser et al. (1994): = (a 2 — l)/(a 2 — a 2 r ) 

where a = r 0 /r max . 

Another crucial point is relevant to the application of 
weak lensing inversion technique to the central region of clus¬ 
ter close to the Einstein radius, where the actual compar¬ 
ison of strong lensing and weak lensing determined cluster 
masses becomes possible. However, because the giant arc 
traces approximately the Einstein radius and because the 
weak lensing method can be marginally applicable to such 
a small radius, large uncertainties could be introduced into 
the shear measurements within the positions of arcs/arclets 
(e.g. Seitz & Schneider 1995). This even does not account 
for the possible contamination of central cluster galaxies due 
to the small number of background faint galaxies near the 
Einstein radius. From Fig.2 it is easily recognized that the 
innermost radii of the weak lensing measurements are al¬ 
ways larger than or equal to the arc radii. Therefore, one 
cannot exclude the possibility that the disagreement of the 
strong and weak inferred cluster masses arises simply from 
the inappropriate application of the week lensing inversion 
technique to the central regimes of clusters, although it has 
been shown that the cluster mass is only slightly underes¬ 
timated at the cluster center if a general method for the 
nonlinear cluster mass reconstruction is used (Kaiser 1995; 
Schneider & Seitz 1995; Seitz & Schneider 1995). 

4.3 X-ray analysis: underestimates cluster mass ? 

An accurate estimate of cluster mass from the X-ray mea¬ 
surement depends on our understanding of the gas distribu¬ 
tion and its dynamical state in cluster, which is closely con¬ 
nected to how significant the X-ray emitting gas deviates 
from the state of hydrostatic equipartition with the grav¬ 
itational potential of the whole cluster. This argument is 
twofold: (l)The complex dynamical activities in cluster may 
only take place locally and on small scales such as the sub¬ 
structure merging and irregular temperature patterns, while 
the cluster as a whole can be regarded as a violently relaxed 
system; (2)Cluster is still in forming stage and thereby, can¬ 
not be modeled as a virialized system at all. 

I 11 the first circumstance, the Jeans equation may be 
safely applied to large radius but fails at small scales. Such 
a dynamical model can account for the excellent agreement 
between weak lensing and X-ray cluster masses at relatively 
large radius, and the mass discrepancy between X-ray and 
strong lensing measurements within the central cores of clus¬ 
ters. Further evidences for this scenario come from the over¬ 
all a ga i-T relationship for the lensing clusters (Fig.l) and 
the recent observations of cluster abundances and other dy¬ 
namical properties at redshift out to z ~ 0.8 (Mushotzky, 
& Scharf 1997; Bahcall et al. 1997; Rosati et al. 1998; etc.), 
i.e., the significant cosmic and dynamical evolutions do not 
play a dominant role for clusters of galaxies as a whole since 
2 ~ 0.8. There are a number of mechanisms that could affect 
the determinations of cluster masses on small scales based 


on the X-ray observations, among which the presence or lack 
of the cooling flows inside the cores of clusters is likely to be 
the major source of uncertainties (Allen 1998). If the core 
radius of the X-ray emitting gas profile is overestimated due 
to the contamination of these central dynamical activities, 
we can indeed reconcile the strong lensing determined clus¬ 
ter masses with the X-ray measurements by substantially 
reducing the X-ray core radius (Fig.5). Recent analyses by 
Markevitch (1997) and Girardi et al. (1997b) for A1689 and 
A2218, based on the high resolution X-ray/optical images, 
have revealed that the local dynamical activities like the on¬ 
going subcluster mergers in the core regions probably make 
the hydrostatic mass estimate inapplicable. Additional sup¬ 
ports come from the numerical simulations of cluster forma¬ 
tion and evolution. For instances, Bartelmann & Steinmetz 
(1996) showed that utilizing a (3 model for a dynamical ac¬ 
tive cluster may lead to an average underestimate of the true 
cluster mass by a factor of ~ 40%, and strong lensing pref¬ 
erentially selects clusters that are dynamically more active 
than the average. From an extensive analysis of the (3 model 
as a cluster mass estimator, Evrard, Metzler & Navarro 
(1996) concluded that while the mass estimates based on 
the hydrostatic, isothermal (3 model are remarkably accu¬ 
rate, the ratio of the estimated cluster mass in terms of the 
(3 model to the true cluster mass increases with increasing 
radius, with the true cluster mass being underestimated by 
the (3 model toward cluster center when flo = 1, and merging 
rates are increased at low redshift. This is fairly consistent 
with our speculation ( 1 ). 

The second possibility is of interest when a large 
universe (e.g. = 1 ) is preferred, in which clusters of 

galaxies formed very recently and may still be in a stage of 
dramatic dynamical evolution. This may essentially invali¬ 
date the application of the Jeans equation without consid¬ 
ering the infall motions of the materials inside and around 
the clusters. The coincidence between the weak lensing and 
X-ray determined cluster masses (see Fig.4) does not ensure 
that both methods would not result in an underestimate of 
total cluster masses. It is somewhat unfortunate if we have 
to accept this scenario, though there have been no convinc¬ 
ing observations so far to confirm this argument. 

4.4 Velocity dispersion of galaxies: a good 
indicator of cluster ? 

Unlike the diffuse X-ray gas, cluster galaxies are unaffected 
by the presence of the (non)cooling flows and the nonthermal 
pressure such as magnetic field and therefore, are probably 
a better indicator of mass distributions of clusters. How¬ 
ever, a precise measurement of the galaxy velocity dispersion 
to large radius from the cluster center is difficult. Employ¬ 
ment of the Jeans equation for galaxies among the ensem¬ 
ble clusters including lensing ones at intermediate redshifts 
have become possible only recently (Carlberg, Yee, & Elling- 
son 1997). On the other hand, the virial and lensing cluster 
masses are basically consistent with each other (see Fig.7), 
though the two methods probe very different regions of clus¬ 
ters. 

Attempts have been made to fit the lensing derived mass 
profile with either a singular or softened isothermal sphere 
model for many individual lensing clusters (e.g. Fahlman 
et al. 1994; Tyson & Fischer 1995; Squires et al. 1997a,b; 
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Luppino & Kaiser 1997; Fischer & Tyson 1998; Tyson, 
Kochanski, & Dell’Antonio 1998; etc), which have neverthe¬ 
less yielded a controversial result. From our statistical study 
(see Fig.6), an uncomfortably small core radius of ~ 0 may 
be required for an isothermal dark matter distribution if the 
strong lensing data alone are used, while the weak lensing re¬ 
sults allow the core size of dark matter profile to be as large 
as the one for luminous matter distributions. A combina¬ 
tion of the strong and weak lensing analyses suggests that 
an isothermal sphere model with a compact but non-zero 
core and characterized by the observed velocity dispersion 
of galaxies a ga i is likely to provide a good description of total 
mass distribution of cluster. This is consistent with the first 
high resolution mass map of CL0024+1654 through param¬ 
eter inversion of the multiple images of a background galaxy 
(Tyson et al. 1998), which has detected the presence of a soft 
core of ~ 0.066 Mpc. Alternatively, if the true mass distri¬ 
bution of cluster is close to an isothermal sphere model, the 
statistical agreement of the lensing measured cluster masses 
with those expected from the isothermal model indicates 
that there is no strong bias between the velocity of the dark 
matter particles and of the galaxies. In other words, veloc¬ 
ity dispersion of galaxies seems to be a good indicator of the 
gravitating mass of cluster. 


The cosmological implications of this work are as fol¬ 
lows: The mass-to-light ratio M/L and the baryon fraction 
fb revealed by the current optical and X-ray observations 
of clusters are indeed reliable indicators for the overall mat¬ 
ter composition of clusters and thereby the universe, though 
these quantities may have rather a large uncertainty on small 
scales of ~ 100 kpc. As a result, we have to accept a low 
mass density universe of Q. m « 0.1-0.4, as suggested by the 
two independent measurements of the cluster matter com¬ 
positions, the ratio M/L [(dark+luminous)/luminous] (e.g. 
Bahcall, Lubin, & Dorman, 1995) and the baryon fraction 
fb [baryon/(baryon+nonbaryon)] (e.g. White et al. 1993). If 
so, we should be aware that the fraction of baryonic matter 
in the universe is not small at all: ~ 20-30% of the matter 
in the universe is luckily visible ! 
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Figure 1. The cr ga i-T relationship of the 20 lensing clusters in Table 1 and Table 2, for which both velocity dispersion <J ga i and 
temperature T are observationally determined. The dashed line is the best-fit to the data: (cr ga i /km s —1 ) = 10 2-57=1=0 13 (fcT/keV) 0,59=,=0,14 . 
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Figure 2. Comparison of the strong and weak lensing determined cluster masses among individual clusters. The open and filled symbols 
represent the strong and weak lensing results, respectively. Note that the giant arc was used for the calibration of the weak lensing 
measurements in A1689, for which we have also displayed the recent result (open asterisk) obtained from the measurement of the deficit 
of red galaxy population behind A1689 for comparison (Taylor et al. 1998). The uncertainties of the strong lensing results are not shown, 
which turn to be quite small when the positions and redshifts of the arcs are observationally determined. 
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Figure 3. The strong and weak lensing measured cluster masses are plotted against the cluster radii for all the data sets in Table 1 and 
Table 2. It appears that the gravitating masses revealed by the arclike images systematically exceed those derived from the inversion of 
the weakly distorted images of background galaxies around clusters. 
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Figure 4. Gravitational lensing determined cluster masses mi ens versus the hydrostatic masses m X ray given by the X-ray diffuse gas. 
Only those clusters in Table 1 and Table 2 whose temperatures are spectroscopically measured are shown. A conventional isothermal 
/3 model with (3fi t = 2/3 and core radius r xray:C = 0.25 Mpc is adopted for the gas distribution. It is remarkable that m xray agree 
essentially with the weak lensing results (a), while an excellent agreement is reached when cluster A2163 is excluded (b). The horizontal 
error bars only reflect the uncertainties in the measurements of temperature. The dashed line is not a fit to the data but assumed that 
TDlens : Dfl xray • 
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Figure 6. Gravitational lensing determined cluster masses mi ens are plotted against the theoretically expected results m op t from a 
softened isothermal sphere model for the total mass distributions of clusters, which is characterized by the velocity dispersion cr ga i of 
galaxies and a core radius r ^ ar *. >c . (a) r^ark^ = 0-25 Mpc and (b) rd ar k,c = 0-025 Mpc. Only the clusters whose <J ga i are observationally 
measured are shown, and the horizontal error bars represent the uncertainties in cr ga i. The dashed line is not a fit to the data but assumed 
that rrii ens = m 0 pt . The dispersions in mi ens can be considerably reduced if the weak lensing data of A2163 and RXJ1716 are excluded. 


© 0000 RAS, MNRAS 000, 000-000 




16 Wu, Chiueh, Fang and Xue 



r (Mpc) 


Figure 7. A comparison of the cluster masses derived from gravitational lensing and virial theorem. At sufficiently large radius, the 
deviation of the projected masses from the spatial (virial) values should become negligible. According to the present data sets, it is 
unlikely that there are significant differences between the two mass estimators. The virial masses are from Carlberg et al. (1996) and 
Girardi et al. (1997a). 
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